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ABSTRACT 


Both numerical and experimental analyses were performed 
to investigate underwater shock propagation and the induced 
nonlinear response of cylindrical shells with end caps. The 
cylinders were subjected to shocks from explosive charges at 
12 inches (near-field) and 28 feet (far-field) from the 
cylinder. An underwater shock test was also performed with the 
far-field explosion. The numerical results were compared with 
the experimental data. Stresses and strains occurring in the 
Structure as well as the pressure in the water were studied. 
The far-field explosion caused the largest circumferential 
deformations close to both end plates and an accordion 
Sscctldatory motion of the cylindrical shell. The near-field 
explosion caused severe plastic deformation in the 
neighborhood of the closest end plate to the charge. The 
stiffeners had, as expected, a larger effect оп the 
circumferential stresses than on the longitudinal stresses. 
The meaSured and calculated strains agreed well qualitatively 


Near the remote end plate from the charge. 
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I. INTRODUCTION 


Because of the U.S. Navy’s high interest in the underwater 
shock hardening effects on surface ships and submarines, this 
research intends to provide more insight into the response of 
a submerged vessel subjected to end-on underwater shock. 

A continuous research has been taken at the Naval 
Peseqracuate - “School го provide. more. insight “into: “the 
deformation and catastrophic failure of surface and subsurface 
hulls. With simple cylindrical shells as a starting point, the 
Study will then be extended to structures with more complex 
material and geometric properties as the methods апа 
predictions improve. 

some of the previous studies in this subject are listed in 
weterences 1-3. These studies have served as building blocks 
for the current research into dynamic response of cylindrical 
S@ells со underwater shock. The mob jective of this study is to 
provide insight into the end-on shock dynamic response of 
Cylindrical shells by using numerical "and experimental 
techniques. 

An unstiffened cylinder subjected to a far-field explosion 
was investigated using both numerical and experimental 
techniques. In addition, both unstiffened and ring-stiffened 


cylinders subjected to a near-field explosion were studied 


numerically. Both the finite element and the boundan, етешеши 


methods were utilized for the numerical study. 


II. ANALYSIS AND EXPERIMENT 


A. NUMERICAL ANALYSIS 

Mie mumetical, Study ef the benavior of cylindrical shel ls 
loaded by underwater explosion was carried out by using the 
finite element and boundary element methods. 

For a problem with a three dimensional domain, the finite 
element method generates a three dimensional discretization of 
the entire domain whereas the boundary element method 
discretizes the surface boundary of the domain with a two 
dimensional grid. The boundary element method reduces 
Significantly the number of elements required to model the 
problem by using a two dimensional mesh. Furthermore, the 
boundary element method surpasses the finite element method in 
computing tractions because these tractions are treated as 
Peimary but not secondary unknowns. Rel. 4]. The matrix 
generated by the boundary element method is generally fully 
populated while the matrix generated by the finite element 
method is usually narrowly banded. 

To study the propagation of the explosive pressure wave 
through the acoustic medium and its subsequent interaction 
with the cylindrical shell, a finite element analysis program 
called VEC/DYNA3D [Ref. 5] was used. This program analyzes the 


deformations of solids which usually are inelastic in nature. 


An alternative for modelling the acoustic medium is by using 
the boundary element method instead of using the finite 
element method. In doing so, the number of elements is greatly 
reduced. The boundary element analysis method program called 
USA (Underwater Shock Analyzer) [Ref. 6], was used to compute 
transient responses of submerged structures to acoustic shock 
waves. 

In USA, the fluid-structure interaction was handled using 
the Doubly Asymptotic Approximation (DAA). The differential 
equation of motion describing the structure response is given 


below. 
MOX + Cea E (1) 


where M,, C,, and K, are the mass, damping and stiffness 
matrices respectively and x, x, xX are the displacement, 
velocity and acceleration vectors of the Struct Ua 
respectively. The excitation for a submerged structure 


subjected to an acoustic wave is given below. 


n c M DEM (2) 


5 а 


where -GA, (P, * P,) is the force vector due to the fluid- 
structure interaction, f 15 the force vector applied о ыш 
dry structure, P. is the incident pressure, P, is the scattered 
pressure, G is the transformation matrix relating fluid and 
structure nodal forces and A, is the diagonal area matrix 


associated with the fluid mesh. 


Dom ru cra cstudl he f$lüsd-cstructure'.snteraction “and 
especially to relate the scattered wave pressure to velocity 
over the wet ее. Фал was InN 7991) "The DAA solution 
approaches the exact solutions for both early time and late 
time responses. The early time response is the high frequency 
response and the late time response is the low frequency 
response. The DAA is given by the following equation. 

рр O. (3) 

Ем. 15 сре Symmetric fluid mass matrix for the wert- 
Surface fluid mesh, u, is the vector of scattered fluid 
різно і сердечно мо вес» nomort' to twecstctrucbure S web surtdgcej vp 
movwene density Of спе fluid, c 15 the sonic speed in tne fluid 
and A. is the symmetric fluid mass matrix for the wet surface 
fluid mesh. This matrix is created using the boundary element 
method. 

The high frequency approximation, i.e. plane wave 
approximation implies that |P,| >>|P,| where P, is the time 
derivative of the acoustic pressure. Therefore, Equation (3) 


шейоопуегіед 1пСсо Спе following equation. 
P, = peu, (4) 


IRE вом rrequency B-apprOximation, те. Virtual, “mass 
approximation implies that |P,| <<|P,| . The concept of virtual 


mass stems from the need to include a mass of fluid 


surrounding that structure at the Mow frequency moon or a 


Structure. Hence, Equation (3) is modified as seen below. 
Ape = ane (5) 


Equation (3) denotes the first order of the Doubly 
Asymptotic Approximations (DAA1). The second order of the 
Doubly Asymptotic Approximations, DAA2, was created to improve 
the solution of DAA1 for intermediate times and to correct for 
curvature of the surface of the structure [Ref. 7]. 

The finite element model of both structure and fluid has 
the advantage of presenting the pressure wave propagation in 
the fluid and the fluid-structure interaction. Оп спе ов Е 
hand, it also requires a large number of elements. Using the 
finite element model for the structure and the boundary 
element model for the fluid reduces the number of degrees Gf 
freedom in the system because the fluid domain is usually much 
larger than the structural dimension. However, it cannot 
present the propagation of the pressure wave across the fluid. 

The post-processing of the VEC/DYNA3D and USA programs was 
done using LS-TAURUS [Ref. 8]. This post-processor generates 
pressure, strain and stress contours superposed on the mesh or 


also generates element and node time history responses. 


B. NUMERICAL MODELLING 
To study the deformation of a cylindrical shell subjected 


to an end-on shock, two numerical models were created. The 


first model was designed to study the propagation of the 
pressure wave from the explosive to the cylinder and its early 
interaction with the cylindrical shell. The second model was 
designed to study stresses and strains of the cylindrical 
shell without analyzing the fluid or exploSive around it. 

To generate the first model, the pre-processor, INGRID 
[Ref. 9], generated the finite element meshes for the 
explosive, fluid and cylindrical shell and VEC/DYNASD computed 
their dynamic response. 

The computational effort was minimized by creating a 
quarter model, possible only because of the symmetric geometry 
of the charge and cylindrical shell locations. Appropriate 
symmetric boundary conditions were applied to the problem. 

The explosive was modeled with a fine mesh of 416 elements 
to avoid non-spherical propagation and severe distortion in 
the finite element mesh as a result of the expansion of the 
explosive in the fluid. 

To model the explosion, the Jones-Wilkins-Lee (JWL) 
equation of state was invoked to describe the pressure-volume- 


eni roy behavior ОТ high, explosives. [RBef. 10] 


P= ali-_”_ b-®” + Bly -_Y bevy WE (6) 
7 ку V 


шнезе А, В, апа Є аге linear coefficients given in units of 





Mbar. R1, R2 and w are nonlinear coefficients, V is the volume 


of detonation products divided by the volume of undetonated 


high explosive, P is the pressure, and E is the detonation 
energy per ünit volumes in (Moark ni E 

The Gruneisen equation of state was used to define the 
pressure for compressed materials, in this case water. This 
equation of state is provided ani ker OE 

The model including the charge, water and cylinder had 
11808 elements which consequently required large amounts of 
storage. Finite element meshes are given in Figure 2.1. 
Therefore, the charge was constrained to be close to the 
target cylinder to minimize the elements between explosive and 
cylinder and still be able to store the shock wave propagation 
and stress wave effect information. The cylindrical shell of 
this finite element model was modified to add ring stiffeners. 
One model had one ring stiffener located halfway between the 
cylinder end plates. The second model had two ring stiffeners 
equidistant between the two end plates. A view of these two 
models is shown in Figure 2.2. The additional stiffeners only 
added an additional 16 elements per stiffener to the total. 
The ring stiffeners were added to study their effect in stress 
wave propagation and reinforcement properties. 

The other model used in this research was designed to 
Study the stresses and strains of the cylinder for comparison 
with experimental results of an actual underwater shock test. 
The model dynamic response was analyzed using both VEC/DYNA3D 
and USA. The finite element method was used for the 


cylindrical shell and the boundary element method was used for 
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(d) Quarter Fluid Mesh 


Figure 2.1. Finite Element Meshes for Near-Field Explosion 
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(a) One-Stiffener Model 
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(b) Two-Stiffener Model 


Figure 2.2. Stiffened Cylinders with End Plates 
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BHEENtluid-etructure interaction ~ This model “consisted. oF а 
cylindrical shell capped at both ends by two end plates. The 
eyiaimarical shell. had 512 elements and the endplates had 48 
elements each. The full model had 608 total elements. A 
diagram of this model is presented іп Figure’ 2.3. The 
different size of the cylindrical shell elements was generated 
to get a better solution in the strain gage locations of the 
experiment. These locations will be presented when the 
experimental procedure is explained. 

A study of grid independence and time step instability was 
performed to verify that the choice of the element sizes and 
time steps chosen were not affecting the solution. To analyze 
grid independence, three different discretizations were made 
over a four inch segment of the cylinder. The discretization 
was made in the longitudinal direction due to the Symmetry of 
the problem. The coarse model had two elements of two inches 
each in the axial direction. The finer model had six elements 
of 0.66 inches each in the same direction. The finest mesh 
model had eighteen elements of 0.22 inches each. The three 
discretizations are shown in Figure 2.4. Comparison of the 
three models revealed considerable differences between the 
coarse model and the two finer mesh models. The six and 18 
element models had almost identical results. As a result of 
this comparison, the six element model was chosen as a good 
estimate of the size of the elements to use for the actual 


model. 
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(a) End plate 
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(b) Cylindrical shell 


Figure 2.3. Finite element meshes of a cylinder subject to a 
far-field explosion 
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(a) Coarse model (2 inch elements) 
(2 elements in 4 inch span) 
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(b) Finer model (0.66 inch elements) 
(6 elements in 4 inch span) 
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(C) Finest model (0.22 inch elements) 
(18 elements in 4 inch span) 


Figure 2.4. Finite element discretization over 4 inch span at 
one end of cylinder 
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Another area of concern was the effect the time step woul 
have on the solution. An appropriate method to choose a stable 
time step size was the Courant-Friederichs-Lewy Criterion. The 


equation for this criterion is givent o O 


۸ ک٥‎ ۴ (7) 


where At is the time step size, @ is the length of the 
element, and c is the sonic speed in 6061-T6 Aluminum. 

Two different time step Sizes =weremcnocen, کت‎ 0e NE لل‎ 
2x10e^! seconds. These time steps satisfied the above mentioned 
criterion. The resulting solutions for the two time steps were 
almost identical. Therefore, 4Ах10е’ зесопа$ миаз свозейп асы 
time step for this problem since this value was between the 


two chosen or tke الال‎ ۹ 1 Ес 


С. UNDERWATER EXPLOSION TEST 

The underwater explosion test was performed at Dynamic 
Testing Inc. (DTI) facilities in RBustburg, Viren qub. 
facility had a quarry that had been filled with water for use 
in underwater shock tests. The water depth was approximately 
130 feet at the location of the test which was deep enough to 
allow the study of the dynamic response of the cylinder prior 
to the arrival of the reflected shock wave from the bottom. 

The cylinder was placed 12 feet below the water's surface 
and held in place by a crane with pendants attached to both 


front and rear end plates. The rig attachment can be observed 
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томе оне 60 “pound -charce ог НВХІі was also plaeed 1Z 
feet below the surface and aligned with the cylinder with a 
Sear wire from"the charge float to the crane rig. Figure 2:6. 
gives a good view of the arrangement. The 12 foot depth was 
chosen so that the bubble generated by the explosion would 
vent to the surface prior to encountering the cylinder. The 
test geometry is shown in Figure 2.7. The explosive charge was 
ati ated ру а radio device and the plume is pictured in 
тшше 2,6. 

The strain gages used for this test were of type CEA-06- 
250UW-350. These are general purpose strain gages with an 
optimum operating range of + 1500 micro strain and are used 
for both static and dynamic test measurements. They were 
attached to the cylinder using a M bond 200. There were a 
total of seven strain gages placed at locations A, В апа С, as 
seen in Figure 2.9, per axis for a total of 14 strain gages. 
The cylinder was oriented so that the gages at C would be 
closest to the explosive charge. A pressure probe was also 
placed 28 feet from the cylinder to measure the free field 
pressure. 

The cylinder used for the underwater shock test was 
па зассиагеа from 6061-T6 Aluminum. This alloy is primarily an 
Aluminum-Magnesium-Silicon alloy. The TG denomination 
indicates it was solution treated and artificially aged [Ref. 
ши ee Live’) cylinder consrsted.-oroas*Ndquarter INCI “thick 


ey lindrical shell and two one-inch thyck circular end plates 


1 


رب 


ре 
мі 
255495 





igging 


Cylinder and Crane R 


Figure 2.5. 
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Figure 2.6. Test Geometry at the Surface 
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Figure 2.9. Cylinder Geometry 
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ой іп күсте 229: Тһе welding of the two end plates to 
the cylindrical shell was done at the Naval Postgraduate 
School facilities. There was a concern in the welding process 
of the two circular end plates to the shell because of the 
heat generated at the weld and its effect on the alloy 
morphology close to the weld. However, the strain readings 
were taken far enough from the heat affected zone to have any 
measurable effect on the readings. The welding was done using 
tungsten inert gas (TIG). This procedure has been recommended 
by most expert welders [Ref. 12]. 

The cylinder weighed 60.5 pounds and tensile tests were 
done to verify that the material properties were close to the 
nominal properties of 6061-T6 Aluminum. This test can be seen 
in the Appendix. The value of Young’s Modulus was 10800 ksi 
and the yield strength was approximately 43 ksi. 

Post-shot observation showed no visible deformations on 
the cylinder as seen in Figure 2.10. All strain gages were 
well fixed when uncovering the bonding material. Only strain 
gage Cl in Figure 2.9. was found wet when uncovered suggesting 


the possibility of water insertion. 
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III. RESULTS AND DISCUSSIONS 


A. NUMERICAL RESULTS OF FAR-FIELD EXPLOSION 

A cylinder was subjected to an end-on explosion with a 
Standoff distance. The sketch of the configuration was 
presented in previous Figure 2.6. The deformations of the 
cylinder are shown in Figure 3.1 at different times as the 
Shock wave propagates along the axial direction of the 
cylinder. Because of the symmetric loading, the deformation 
was axially symmetric. In order to visualize the deformations 
more clearly, the actual deformations were magnified by the 
scale factor shown in Figure 3.1. 

As the shock pressure wave hit one end plate, which is at 
che nearest location from the charge, the compressive stress 
wave propagated from the end plate to the other end plate 
through the aluminum alloy cylinder with a faster speed of 
sound than the shock pressure wave propagating through the 
Surrounding water medium. The nearer end plate will be called 
the near plate while the other end plate will be called the 
remote end plate in the following discussion. The shock load 
caused initially a localized circumferential deformation close 
to the near end plate as shown in Figure 3.lc. This localized 
deformation remained as a permanent deformation. As the shock 


pressure wave propagated toward the remote end plate, it 
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گ تک دگ penes‏ 










ІШІ‏ هف 
а ША ААА ИА‏ 


ар 1 БЕ О ۱ ره د‎ у 


mc ан унцій и ————M——ÉEÓÀ 


(d) time = .15960E-03 sec 










ش 2 تت FEET‏ 


ЕНЕ 
СОД 111, 


tt بل‎ ppt 








(e) time = .23960E-03 sec 





٢ 28 27:2 22 به مه مه‎ ١ ممه‎ ТГ هم‎ ERR RUN ro 


Е ЕЕЕЕЕЕЕ 0009 
"ИТТЕЙ 


L| 
ШІДЕН c HHH 
АЕ) 


На ER) 








(Ё) time = .47960E-03 sec 


disp. scale factor = 1.000E+01 


(cont). Deformation Time Histories 


25 


ART ore Чү Ола ше 
ОЮН 


СШ EEEO | 
TT 


аана у) 

















А 





51960Е-03 sec 












(Cee =: 
TT Пет E US 
ГЕ 8 a 






ашу 


ptc ле КА сили ыы л отте -1- +---<-------х4--2.л-жх-1х44 47“ у--- 






(1) time = .11596E-02 sec 


disp. scale factor = 1.000E+1 


Figure 3.1 (cont.) Deformation Time Histories 


26 


ШІСНСЕІЕП a ~radial contraction «оғ tbe cylinder wall апа 
subsequently a recovery of it. After the shock pressure wave 
passed the remote end plate, there was also a localized 
Circumferential deformation near the remote end plate as shown 
in Figure 3.1h. The overall steady state deformation of the 
cylinder was nearly symmetric about the center plane which is 
located at an equal distance between the two end plates. At 
this stage, the cylinder had the deformation near the two end 
plates, and the rest of the cylinder had very little radial 
deformation. The axial deformation of the cylinder had an 
doceo bn mode as expected. Figure 3.2 1$ the plot of. axial 
velocities at the centers of the two end plates. The two 
velocities had a phase difference of 180 degrees. The phase 
difference indicated the accordion mode. 

The initial localized circumferential deformation near the 
end plates may be explained as follows. The axisymmetric 


deformation of a shell has the following governing equation. 





4 2 
peru жан Е р (8) 
Gs ах“ а: 


where w is the radial deflection, x is the axial direction, D 
Ше (пе flexural rigidity of the shell, Е 15 the elastic 
modulus, h is the wall thickness of the shell, a is the shell 
Pacilus, and Е ап й p are the axial тоза and the radial load, 
respectively. Therefore, the effect of the inertia term on the 
deformation was negligible. The axial load applied on the 


cylinder due to the shock pressure was lower than the critical 
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ОиСТ Ба dosc of that cylinder trom the linear buckling 
theory. In this case, the deformation of the cylinder was of 
a decaying exponential form and it damped out rapidly if the 
cylinder was not short [Ref. 13]. The very stiff end plates 
suppressed the deformation very near the end plates. This 
resulted in a shift of the locations of the circumferential 
deformation somewhat away from the end plates. The inertial 
force was not included in the equation. When the density of 
the cylinder was varied from its nominal value to a value of 
one order of less magnitude, the same kind of initial 
localized deformation was observed. However, the low density 
caused less permanent localized deformation. 

Two normal strain components, Deu hoop (or 
circumferential) and axial (or longitudinal) strains, were 
computed at some selected locations. The selected points were 
located at the locally deformed zones as well as at the center 
of the cylinder. The locations where the strains were computed 
are illustrated in Figure 3.3. Strain gages were also attached 
to the same locations on the cylinder in the experimental 
study for the comparison with the numerical study. All the 
strains were computed and measured at the outside surface of 
the cylinder. The comparison between the numerical and 
experimental results will be provided in the next section. As 
shown in the strain plots of Figure 3.4, the circumferential 
Strains at locations A and C of Figure 3.3 were very close to 


each other. The strains had initially large compressive peaks 
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and returned to the steady state values which were about half 
the initial peak strain values. Location A had a small 
Circumferential tensile strain before arrival of the shock 
pressure wave. The small tensile strain was caused by the 
longitudinal compressive stress wave which arrived at the 
location through the cylinder. The longitudinal compressive 
Strain caused by the longitudinal stress wave resulted in the 
Circumferential tensile strain because of Poisson’s effect. 
After the shock pressure wave arrived at the location, the 
Circumferential strain became compressive. The circumferential 
Strain at the middle of the cylinder, i.e, location TE s 
Figure 3.3, was much less that those at locations A and C as 
shown in Figure 3.5 because the localized deformations at 
locations A and C induced the larger circumferential strains. 

The initial peak of longitudinal strain was compressive at 
all locations due to the propagation of the longitudimam 
compressive stress wave after the shock pressure hit the near 
end plate. The early arrival of the shock pressure wave to 
location C in Figure 3.3, caused the immediate jump or Tens 
longitudinal strain Ineo tension after the 11111 1 
compression. On the other hand, the subsequent compressive 
stress wave resulted in more compression at locations A and B 
than at location C until the shock pressure wave arrived at 
the locations. At later times the longitudinal strain were 
quite different at locations A and C. The longitudinal stra 


at location A stayed іп tension just after the initial 
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compressive peak while the longitudinal strain ас Iik ati TE 
alternated between tension and compression. However, the 
longitudinal strain was quite sensitive with the locations 
near A and C because of the localized deformation. 
Longitudinal strains at various locations near position C were 
plotted in Figure 3.6. The approximate frequency of the 
longitudinal strains after the initial compression was about 
750 HZ at a11 ه109‎ ене 

The present study showed that the end plates played an 
important role in the deformation of the cylinder сеща 
parametric study was performed to find the effect of the end 
plates on the dynamic response of the cylindrical shel ال‎ 
first study was to find out the effect of the inertial force 
of the end plates. Therefore, the density of the end plates 
was reduced tenfold with the stiffness of the plates remaining 
the same. The light end plates had a larger effect on the 
longitudinal strains than on the circumferential strains. The 
frequency of the longitudinal Strains was approximately 840 Hz 
for the light end plates. This frequency was higher than that 
for the heavy plates as expected. The longitudinal strain 
damped more quickly with the light end plates than with heavy 
end plates as shown in Figure 3.7. The light end plates 
resulted in an approximately 30% decrease one the 
circumferential strain at location C while its effect was 
negligible on the circumferential strain at location A. 


Moreover, the circumferential strain at location B increased 
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about 20 percent for the light end plates as shown in Figure 
DOLO 

The next parametric study was performed to investigate the 
effect of stiffness of the end plates on the deformation к 
the cylinder. The elastic modulus of the plate was reduced 
tenfold without change of the rest of the material data. The 
less stiff end plates made an effect on not only longitudinal 
Strain but also circumferential strain as shown in Figures 3.9 
and 3.10. The less stiff support of the cylinder at the end 
plates caused more fluctuation in the circumferential strain 
as shown in Figure 3.9. The wave pattern in the longitudinal 
Strain-time history plot was also severely altered by the less 
stiff end plates as shown in Figure 3.10. The longitudinal 
velocities of the centers of the less stiff end plates were 
plotted in Figure 3.11 and showed a very different pattern 
than the accordion oscillation of the stiff plates in ۳٣ ٢ 
3.2. The density change in the previous parametric study did 


no show a different velocity pattern of the end plates. 


B. COMPARISON BETWEEN NUMERICAL AND EXPERIMENTAL RESULTS 

An underwater explosion test was performed at Dynamic 
Testing Inc. facilities in Rustburg, Virginia as described а 
Chapter II. Unfortunately, a pressure gage, placed to measure 
a free-field pressure at the same stand-off distance as the 
cylinder, failed. Therefore, no information was available for 


the shock pressure due to the charge. The previous numerical 
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Simulation was carried out using the empirical equation for 
pressure with the nominal weight of the charge. As a result, 
the following comparison was made in a qualitative sense. In 
addition, several strain gages failed before and during the 
test. All the strain gages to measure circumferential strains 
failed. Only the axial strains were compared here. 

The speed of sound in the water was computed based on the 
arrival time of the pressure wave to the strain gage from the 
explosion and the stand-off distance. The computed speed of 
sound was close to a nominal value of 5000 ft/sec. Then 
locations of the strain gages were shown in Figure 3.3. The 
longitudinal strain at gage location A was compared in Figure 
3.12. All experimental data was filtered out at 2000 Hz low 
pass. Therefore, there was no peak strain with a higher 
frequency than 2000 Hz in the experimental data. Both 
numerical and experimental solutions indicated the initial 
compression and tension at later times. They agreed well 
qualitatively even if there was a mismatch in magnitude. The 
longitudinal strain was quite sensitive to the location. The 
longitudinal strain computed at just one inch away from the 
gage location was compared with the experimental measurement 
in Figure 3.13. Two strain gages, both of which were located 
at location A but were separated in the circumferential 
direction with an angle of 180 degrees, gave a similar 
longitudinal measurement as shown in Figure 3.12. This 


measurement indicated the axisymmetric nature of the 
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deformation. However, the two strain gages, which were placed 
at location C and were separated 180 degrees in the 
circumference, provided quite different longitudinal strains. 
This discrepancy in readings could have been caused due to 
observed water insertion in the strain gage at location Cl. 
Accordingly, the numerical result was not in good 
agreement with the experimental data as shown in Figure 3.14. 
The sensitivity of the location were the strains were computed 
was also evident for location C as shown in Figure 3.15. The 
experiment showed a compressed longitudinal strain while the 
pumerical solution showed an oscillation of the strain from 
compression to tension. A small distance away from location C 
gave a different solution. However, the experimental data was 
Still different from the numerical result. The longitudinal 
Strain at gage location B was compared in Figure 3.16. The 
Bolles wes similar. to that at gage location C. It was пої 
clear at this time what the major cause of the discrepancy 
was. However, this was the first test of a series of 
underwater experiments to be performed. The discrepancy will 
be investigated from the following experiments with more 


available experimental data and numerical simulations. 


C. NUMERICAL RESULTS OF NEAR-FIELD EXPLOSION 
The study of the dynamic response of cylindrical shells 
subjected to a close-in explosion was performed. The stand-off 


distance was 12 inches from one end plate. The overall 
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4.5 Inches from Closest End Plate to 
Explosive 
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Figure 3.14. (cont.) Measured and Computed Longitudinal 
Strains at 4.5 Inches from Closest End Plate 
to Explosive 
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20.75 Inches from Closest End Plate to 
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geometry with finite element meshes is shown in Figure 2.1. A 

two pound explosive charge of pentolite was used in the 

numerical study. The dimension of the cylinder was the same as 

for the far-field explosion except for a length of 43 inches 

instead of 42 inches long and the addition of ring stiffeners. 
l. Pressure Wave Study 

The explosive charge was refined in the numerical 
model to avoid discontinuities of the sphere as it expanded 
with time. The pressure wave propagated spherically through 
the finite element mesh as shown in Figure 3.17. 

The wave propagation was studied especially at the 
first interaction between the pressure wave and one end plate. 
The cylinder was accelerated when the wave first interacted 
with the end plate. At about 0.3 milliseconds the fluid and 
end plate velocities break away for a short period. The fluid 
velocity was lower than that for the shell. Consequently, 
tension would be induced in the fluid but since fluids cannot 
experience tension, the water particles break away creating a 
vacuum. This effect, known as hull cavitation, can be observed 
Mee Toure 3.408. 

2. Ring Stiffener Study 

The three cylinders studied were an unstiffened 
Cylinder, a one-ring stiffener cylinder, апа а two-ring 
stiffened cylinder. The ring stiffeners were located at equal 


intervals between the end plates. In other words, one ring 
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Em 


stiffener was placed at the center of the cylinder as seen 
previously in Figure 2.2a and the two ring stiffeners were 
placed at one third distances from the end plates respectively 
as seen previously in Figure 2.2b. The stiffeners were one 
inch high and 0.25 inches thick. The neminal mabertal prop rS] 
data of a 6061-T6 Aluminum alloy was used for this study. 

The ring stiffened cylinders were compared to the 
unstiffened cylinder to provide insight into the effectiveness 
of the ring stiffeners in affecting stress wave propagation 
and deformation of cylinders. The numerical model was labeled 
in specific locations where the dynamic response was of 
interest. These locations are indicated in Figure 3.19 and 
they include the largely deformed area at the front end of the 
cylinder, F, the stiffener locati@ms: Sl, S2 17 53 ane ве: 
location close to the remote end plate, R. 

The initial shock pressure produced a severe local 
deformation about three inches from the front end plate for 
all three cases of cylinders as shown in Figure 3.20. The 
cylinder yielded identically in all cases because the 
stiffeners had no effect on the initial deformations near the 
front end plate until the pressure wave passed through the 
stiffeners. Comparison of the circumferential strains at 
location S1 between the unstiffened and two-stiffener models 
was given in Figure 3.21. In the figure legends A and B 
denoted shell elements located just before and after the 


stiffener. The stiffener reduced the compressive peak strain 
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FE: . 2.6 inches trom end-piate closest to charge 
Sl: First stiffener of two-stiffener model 

S2: Stiffener of one-stiffener model 

53: Second stiffener of two-stiffener model 


R: 40.4 inches from end-plate closest to charge 
Figure 3.19. Near-Field Model Locations Chosen for 
Analysis 


at the location of the stiffener as expected. However, the 
unstiffened cylinder had a larger relief of the compressive 
strain than the two-stiffened cylinder. Comparing the same at 
location S3 indicated that the stiffener not only reduced the 
compresSive peak strain, but also relieved the strain more at 
later times as shown in Figure 3.22. The circumferential 
strains at the two stiffeners were plotted in Figure 3.23. The 
front stiffener had a much larger strain than the back 
stiffener. The large plastic strain at the front stiffener 
restrained the recovery of the circumferential strain at the 
adjacent shell elements at later times compared to the 
unstiffened cylinder as shown in Figure 3.22. 

Comparison of the unstiffened and one-stiffener 
cylinders revealed that the stiffener caused a large relief of 
the circumferential strain. The center stiffener also caused 
a big difference in the circumferential strains between the 
shell elements just before and after the stiffener. The strain 
was greatly reduced after the stiffener as seen in Figure 
3.24. The longitudinal strains were compared in Figures 3.25 
and 3.26 between two-stiffened and unstiffened. Figure 3.25 
was the comparison at location S1 and Figure 3.26 at location 
S3. The ring stiffener at location S1 altered the longitudinal 
Strain before and after the stiffener as shown in Figure 3.25. 
The stiffener at location S1 induced a higher longitudinal 
strain at the shell element located just before the stiffener 


compared to that at the same location of the unstiffened 


20 





(b) One-stiffener model 
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(c) Two-stiffener model 


یک مو ووو Lime o‏ 


Figure 3.20. Local Deformations for Near-Field Models 
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(D) One-stiffener model 


Figure 3.21. Circumferential Strains at Both Sides of 
Stiffener Closest to Explosive Charge 
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(b) One-stiffener model 


Figure 3.22. Circumferential Strains at Both Sides of 
Stiffener Farthest from Explosive Charge 
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(a) Stiffener 14.5 in. from end close to explosive 
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(D) Stiffener 28.5 in. from end close to explosive 


Figure 3.23. Circumferential Strains at Both Stiffeners of 
Two-Stiffener Model 
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cylinder. In addition, the stiffener relieved the strain near 
to zero at the shell element located just after the stiffener 
as time elapsed. The effect of the stiffener at location 53 оп 
the longitudinal strain was much smaller compared to that at 
location S1. No accordion mode was observed for all three 


cases of cylinders. 
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(a) Unstiffened model 
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(D) One-stiffener model 


Figure 3.24. Circumferential Strain Comparison for 
Unstiffened and One-Stiffener Models 
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(a) Unstiffened model 
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(b) Two-stitfener model 


Bigure 3.25, Longitudinal Strains (Unstiffened and Two- 


Stiffened) at Stiffener Location Closest to 
Charge 
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(a) Unstiffened model 
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(b) Two-stiffener model 


Figure 3.26. Longitudinal Strains (Unstiffened and Two- 
Stiffened) at Stiffener Location Far from 
Charge 
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IV. CONCLUSIONS 


The numerical and experimental study was performed to 
investigate the nonlinear dynamic response of cylindrical 
shells to end-on explosions. Unstiffened and ring-stiffened 
cylinders capped at both ends were subjected to near-field and 
far-field explosions. Both finite element and boundary element 
methods were used for the numerical study. The experiment was 
0111 ه1‎ uscmg' a Oo06l-To aluminum alloy cylinder and a “sixty 
pound HBxX!Icharge with 28 feet stand-off distance. 

The far-field explosion resulted in a nearly Symmetric 
deformed shape of the cylinder about the center plane between 
the two end plates, while the near-field explosion caused very 
unsymmetric deformation about the center plane. Both far-field 
and near-field explosions induced localized deformations near 
the end plates. The localized deformations were located close 
Бой сеп end plates for the far-field explosion and a very 
severe, localized deformation occurred close to the closest 
end plate to the explosive charge for the near-field explosion 
regardless of the existence of ring stiffeners. The accordion 
mode was observed for the cylinder subjected to the far-field 
С шоп Dut not forthe cylinders subjected Со пеаг-Ітеїа 
explosion: 

A parametric study was undertaken to examine the effect of 


end plates on the deformation of the cylinder subject to the 


65 


far-field explosion. A variation of the density and stiffness 
of the end plates caused a significant change in the stress 
wave propagation and the deformation of the cylinder. However, 
the change of stiffness had more significant effects than the 
change of density. 

The circumferential strains under the far-field explos ٢ل‎ 
had large initial compressive peaks and returned to the steady 
state values quickly. The circumferential strains were larger 
near the end plates than around the center due to the local 
deformation. The localized deformations were caused by the 
compressive shock pressure applied on the end plates. The 
longitudinal strains under the far-field explosion were 
compressive at the shell near the end plate closer to the 
charge but tensile at the opposite location. 

The failure of pressure and several strain gages prevented 
a quantitative comparison between the numerical and 
experimental results. A qualitative comparison was possible 
between the two solutions only at a few locations. The 
comparison was better near the remote end plate than the near 
end plate to the charge. A good explanation ҒО “ӘРЕ 
discrepancy was not possible due to the lack of information in 
the experiment. However, this was the first test among a 
series of experiments to Бе performed. More detailed 
information from the next tests may provide a better 


understanding of the discrepancy. 
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The effect of ring stiffeners was larger оп the 
s" cmrerenbpra strain than on they longitudinal strain. The 
ring stiffener located close to the severe local deformation 
lad а лоте pronounced effect on the strains. The ring 
stiffeners in general reduced the circumferential strain 
compared to the unstiffened case and stiffeners also gave more 
recovery of the strain. However, when the stiffeners had a 
large plastic deformation, they allowed less relief of the 


emn Un erential strain: 
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APPENDIX: UNIAXIAL TENSION TEST DATA FOR 6061-T6 ALUMINUM 


A uniaxial tension test was performed on a test specimen 
of 6061-T6 Aluminum to find its material properties. The 
specimen was cut to ASTM E8-69 specifications for rectangular 
tension test specimens. 

A Material Test System (MTS) model 810 was used for the 
test and the gage length was set at two inches for the tensile 
test. The tensile test was performed at a rate of 400 sec/inch 
and a lad versus displacement curve was plotted as shown in 
Figure A.1. The MTS model did not have data recording 
capabilities and therefore, some points were taken from the 
graph after the test to plot the appropriate stress-strain 
curve. Some sample data points are listed in Table A.1. 

From the stress-strain curve in Figure A.2, the yrii 
stress was found to be close to 43 kSi uSing a 0.2 percent 


offset and the modulus of elasticity was 10800 ksi. 
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Figure A.1. Load-Displacement Curve for 6061-T6 Aluminum 
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